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The thermal conductivity k, of the iron-arsenide superconductor Ba(Fei_ I Co 2 ,)2As2 was measured 
down to 50 mK for a heat current parallel (k c ) and perpendicular (/c a ) to the tetragonal c axis, for 
seven Co concentrations from underdoped to overdoped regions of the phase diagram (0.038 < x < 
0.127). A residual linear term k c o/T is observed in the T — > limit when the current is along 
the c axis, revealing the presence of nodes in the gap. Because the nodes appear as x moves away 
from the concentration of maximal T c , they must be accidental, not imposed by symmetry, and 
are therefore compatible with an s± state, for example. The fact that the in-plane residual linear 
term K a o/T is negligible at all x implies that the nodes are located in regions of the Fermi surface 
that contribute strongly to c-axis conduction and very little to in-plane conduction. Application of 
a moderate magnetic field (e.g. H C 2/4) excites quasiparticles that conduct heat along the a axis 
just as well as the nodal quasiparticles conduct along the c axis. This shows that the gap must be 
very small (but non-zero) in regions of the Fermi surface which contribute significantly to in-plane 
conduction. These findings can be understood in terms of a strong k dependence of the gap A(k) 
which produces nodes on a Fermi surface sheet with pronounced c-axis dispersion and deep minima 
on the remaining, quasi-two-dimensional sheets. 

PACS numbers: 74.25.Fy, 74.20.Rp,74.70.Dd 



I. INTRODUCTION 

The discovery of superconductivity in iron arsenidespJ 
with transition temperatures exceeding 50 Kpl breaks 
the monopoly of cuprates as the only family of high- 
temperature superconductors, and revives the question 
of the pairing mechanism. Because the mechanism is 
intimately related to the symmetry of the order param- 
eter, which is in turn related to the k dependence of the 
gap function A(k), it is important to determine the gap 
structure in the iron-based superconductors, just as it 
was crucial to establish the d-wave symmetry of the gap 
in cuprate superconductors. The gap structure of iron- 
based superconductors has been the subject of numerous 
studies (for recent reviews, see Refs. 1314]) . Here, we focus 
on the material BaFe2As2, in which superconductivity 
can be induced either by applying pressure^ or by various 
chemical substitutions, such as K for Ba (K-Bal22jP or 
Co for Fe (Co-Bal22) In the case of Co-Bal22, single 
crystals have been grown with compositions that cover 
the entire superconducting phase (see Fig. pEQl 

Two sets of experiments on doped BaFe2As2 appear to 
give contradictory information. On the one hand, angle- 
resolved photoemission spectroscopy (ARPES) detects a 
nodeless, isotropic superconducting gap on all sheets of 
the Fermi surface in K-Bal22^S and in optimally-doped 
Co-Bal22p21and tunneling studies in K-Bal22 detect two 
full superconducting gapsP^ The magnitude of the gaps 
in ARPES is largest on Fermi surfaces where a density- 
wave gap develops in the parent compounds This is 



taken as evidence for an s± pairing state driven by an- 
tiferromagnetic correlations LMlSHlZl Q n t h e other hand, 
the penetration depth in Co-Bal22j 18 l 19 l the spin-lattice 
relaxation rate in K-Bal22pHl and the in - plan e thermal 
conductivity in K-Bal22^ and Co-Bal22j2223lfor exam- 
ple, are inconsistent with a gap that is large everywhere 
on the Fermi surface. Note, however, that the evidence 
for deep minima in the g ap is particularly clear in the 
overdoped regime poMES a re gime which has not so far 
been probed by either ARPES or tunneling. 

Another possible explanation for the apparent discrep- 
ancy between the two sets of experimental results is a 
different sensitivity to the c-axis component (k z ) of the 
quasiparticle k vector, taking into account the three- 
dimensional (3D) character of the Fermi surfaceP^mSfll 
Nodes along the c axis were suggested theoretically to 
explain the discrepanc y bet ween ARPES, penetration 
depth and NMR studies EU22I A variation of the gap mag- 
nitude as a function of k z was suggested in experimental 
studies of the neutron resonances in optimally-doped Ni- 
Bal22Plt was also invoked to explain the temperature 
dependence of the penetration depth in Co-Bal22^1 and 
its a — c anisotropy in Ni-Bal22p3 Clearly, it has become 
important to resolve the 3D structure of the supercon- 
ducting gap function in doped BaFe2As2- 

Heat transport measured at very low temperatures is 
one of the few directional bulk probes of the gap struc- 
ture. The existence of a finite residual linear term kq/T 
in the thermal conductivity k(T) as T — > is unambigu- 
ous evidence for the presence of nodes in the gapp^lEU 
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FIG. 1: Phase diagram of Ba(Fei_ :E Co; C )2As2 as a function of 
Co concentration x, showing the orthorhombic phase below T s 
(blue), the antiferromagnetic (AF) phase below Tjv (green), 
and the superconducting phase below T c (small black dots), 
as determined from resistivity, magnetization and heat capac- 
ity dataj^ and from x-rajPl and neutrorP^ data. Red circles 
mark the T c value of seven samples measured in the present 
study (the c-axis samples labelled A in Table I), indicating 
the range of concentrations covered. The vertical dashed line 
at x = 0.06 marks the approximate location of the critical 
concentration x s where at T = the system goes from or- 
thorhombic (Ortho) in the underdoped region (to the left) to 
tetragonal (Tetra) in the overdoped region (to the right). At 
T = 0, the AF phase also ends close to x = 0.06P 1 



and thus by measuring n(T) as a function of direction in 
the crystal, o ne can loc ate the position of nodes on the 
Fermi surface d 4 * 35 * 39 * 40 ! Here, we report measurements of 
heat transport in Ba(Fei_ x Coj;)2As2 for a current direc- 
tion both parallel and perpendicular to the c axis of the 
tetragonal (or orthorhombic) crystal structure. Our main 
finding is a sizable residual linear term kq/T for a cur- 
rent along the c axis, and a negligible one for a current 
perpendicular to it. This implies the presence of nodes 
in the gap in regions of the Fermi surface that dominate 
the c-axis conduction and contribute little to in-plane 
conduction. Our study shows that the gap structure of 
Co-Bal22 depends on the 3D character of the Fermi sur- 
face in a way that varies strongly with x. 



II. EXPERIMENTAL 



A. Samples 
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FIG. 2: Temperature dependence of the thermal conductiv- 
ity k, plotted as k/T vs T 2 , measured along the c axis in a 
Co-Bal22 sample with x = 0.127, using the two-probe tech- 
nique, for three values of the applied field: H — 0.0, 0.05 and 
0.10 T. At H = 0, the total thermal resistance of sample plus 
contacts is dominated by the very high thermal resistance of 
the two superconducting tin contacts. At H — 50 mT, tin is 
no longer superconducting and the thermal resistance of the 
two contacts has become negligible compared to the sample 
resistance. In this case, extrapolation of k/T to T=0 gives us 
almost exactly the residual linear term in the sample's ther- 
mal conductivity. Increasing H slightly beyond 50 mT, for 
example to 100 mT, leads to no further change in the data. 
This shows that measurements in H = 50 mT reveal the in- 
trinsic zero-field behavior of the sample. The inset shows the 
arrangement of the tin (Sn) contacts and the silver (Ag) wires 
on the sample. 



concentration, x, roughly 0.7 times the flux load compo- 
sition (or nominal content). We studied seven compo- 
sitions: underdoped, with x = 0.038, 0.042, and 0.048; 
overdoped, with x = 0.074, 0.108, 0.114, and 0.127. In 
this Article, 'underdoped' and 'overdoped' refer to con- 
centrations respectively below and above the critical con- 
centration x s ~ 0.06 at which the system at T = goes 
from orthorhombic (below) to tetragonal (above) P The 
T c value for each composition is shown on the phase di- 
agram in Fig. [I] A total of twelve c-axis and nine a-axis 
samples were studied; their characteristics are listed in 
Tables I and II, respectively. Three of the a-axis samples 
were the subject of a previous study (0.074-B, 0.108-A, 



and 0.114-B)M 



Single crystals of Ba(Fei_ a; Co a ;)2As2 were grown from 
FeAs:CoAs flux, as described elsewhereP The doping 
level in the crystals was determined by wavelength dis- 
persive electron probe microanalysis, which gave a Co 



B. Two-probe transport measurements 

Thermal conductivity was measured in a standard one- 
heater-two-thermometer technique! 4 ^ The magnetic field 
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TABLE I: Properties of the twelve c-axis samples of Co-Bal22 
used in this study, x is the Co concentration measured by 
wavelength dispersive microprobe analysis. The supercon- 
ducting transition temperature T c is the temperature at which 
the resistivity goes to zero. The values of the upper critical 
field H C 2 needed to suppress superconductivity in Co-Bal22 
at T — > are taken from Refs.l8"land l41l (H C 2(T) is defined as 
the end, or 'offset', of the superconducting drop in p(T) vs H). 
The residual resistivity p c o is obtained by a smooth extrapo- 
lation of the p c (T) data to T = 0, as shown in Fig. [3] The 
normal-state residual linear term in the thermal conductivity, 
Kcn/T, is obtained from the Wiedemann- Franz law applied 
to p c o (see text). The zero-field residual linear term, n c o/T, 
is obtained by extrapolating to T = the zero-field thermal 
conductivity k c with a linear fit to K c /T vs T 2 , as shown in 
Fig. [5] Kco/T is also expressed as a fraction of the normal- 
state k c n/T, denoted kq/kn = (kq/T) /(kn /T). 
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H was applied along the [001] or c-axis direction of the 
crystal structure, which is tetragonal for overdoped sam- 
ples and orthorhombic for underdoped samples at low 
temperatures. Data were taken on warming after having 
cooled in a constant field applied above T c to ensure a 
homogeneous field distribution. 

It is conventional to measure electrical and thermal 
resistance in a four-probe configuration to avoid the con- 
tribution of contact resistances. This is what was done 
for data taken with a current in the basal plane (J || a, 
in the notation app ropria te for the tetragonal phase), as 
described elsewhere! 21 * 22 ^ For a current along the c axis 
(J || c), however, the four-probe technique is difficult be- 
cause of the strong tendency of iron-arsenide crystals to 
exfoliation, which makes it difficult to cut samples thick 
enough in the c direction to attach four contacts! 29 * 30 ^ 
Consequently, c-axis transport was measured using a two- 
probe technique, which is valid provided contact resis- 
tances are much smaller than the sample resistance. 

Contacts to the c-axis samples were made using silver 
wires (of 50 pm diameter) , soldered to the top and bot- 
tom surfaces of the sample with ultrapure tin (see inset 
of Fig. [2]). The contact making and properties are de- 
scribed in detail in Ref. [43j In brief, these contacts are 
characterized by a surface area resistivity in the nfi cm 2 
range, which, for a typical sample size, yields a contact 
resistance below 10 pfi. This is negligible compared to 



TABLE II: Properties of the nine a-axis samples of Co-Bal22 
used in this study, x, T c , H C 2 and ko /kn are defined in Table I. 
The residual resistivity p a o is obtained by a smooth extrap- 
olation of the PaiT) data to T = 0, as shown in Fig. [4] The 
normal-state residual linear term in the thermal conductivity, 
Kbn/T, is obtained from the Wiedemann- Franz law applied 
to p a o. The zero-field residual linear term, K a o/T, is obtained 
by extrapolating to T = the zero-field thermal conductiv- 
ity K a with a power-law fit to K a /T (see text), as shown in 
Fig. [7] Note that the magnitude of n a o /T, whether positive 
or negative, is in all cases lower than the uncertainty in the 
extrapolation (see Fig. [9| . 
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a typical sample resistance in the normal state, of the 
order of 10 mf2. 

Because tin is a superconductor, the thermal resistance 
of the contacts at very low temperature is large. We 
therefore have to apply a small magnetic field to sup- 
press the superconductivity of tin and make it a normal 
metal, with the very low electrical and thermal resistance 
mentioned above. A field of 50 mT is sufficient to do this. 
In Fig. [2j we compare data obtained with H = 0.00, 0.05 
and 0.10 T. The effect of switching off the contact resis- 
tance with the field is clear, and once tin has gone normal, 
the data is independent of a further small increase in H . 
We therefore regard the data taken at H = 0.05 T as 
representative of the zero-field state of the sample. 



C. Electrical resistivity 

In Fig. [3j we show the temperature dependence of the 
electrical resistivity p c (T) of our c-axis samples, mea- 
sured in a two-probe configuration. In all samples, the 
resistivity follows qualitatively the temperature depen- 
dence reported previously.^ Data for the a-axis samples 
are shown in Fig. |4j A smooth extrapolation of p(T) to 
T = yields the residual resistivity po listed in Table I 
for c-axis samples and in Table II for a-axis samples. The 
uncertainty associated with the extrapolation of p(T) to 
T = is approximately ± 5 %. Due to the uncertainty 
in measuring the geometric factor, the absolute value of 
the resistivity has an error bar of approximately ± 20 % 
for a-axis samples and a factor of 2 uncertainty for c-axis 
samplespSHUl The higher p values in the underdoped 
regime are due to a reconstruction of the Fermi surface 
in the antiferromagnetic phase. The residual resistivity 
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FIG. 3: Temperature dependence of the c-axis resistivity 
p c (T) for the twelve c-axis crystals of Ba(Fei_ a; Co 2 ,)2As2 
studied here, with Co concentrations x as indicated. Top and 
bottom panels show underdoped and overdoped samples, re- 
spectively. Circles and triangles of the same color correspond 
to two different crystals at the same doping, respectively la- 
belled A and B (see Table I). The lines show how the data is 
extrapolated to T — 0, to determine the value of the residual 
resistivity p c o, given in Table I. 



FIG. 4: Temperature dependence of the a-axis (in- 
plane) resistivity p a (T) for the nine a-axis crystals of 
Ba(Fei_ I Co ;c )2As2 studied here, with Co concentrations x 
as indicated. Top and bottom panels show underdoped and 
overdoped samples, respectively. Circles and triangles of the 
same color correspond to two different crystals at the same 
doping, respectively labelled A and B (see Table II). The lines 
show how the data is extrapolated to T — 0, to determine the 
value of the residual resistivity p a o, given in Table II. 



Po is used to determine the normal-state thermal conduc- 
tivity kn/T in the T = limit via the Wiedemann- Franz 
law, k n /T = L /p a , where L = 2.45 x 1(T 8 W fi / K 2 . 
Because the same contacts are used for electrical and 
thermal measurements, the relative geometric-factor un- 
certainty between the measured n and this electrically- 
determined kn is minimal. 



III. RESULTS 
A. Heat transport in the c direction 

The thermal conductivity of solids is the sum of elec- 
tronic and phononic contributions: k — n e + n p . In the 
T = limit, the electronic conductivity is linear in tem- 



perature: K e oc T . In practice, the way to extract n e is 
to extrapolate k/T to T = 0, and thu s obtain the purely 
electronic residual linear term, kq/tW ^ 2 ^ If one can 
neglect electron-phonon scattering, as one usually can 
deep in the superconducting state, then the mean free 
path of phonons as T — > is controlled by the sample 
boundaries. If those boundaries are rough, the scatter- 
ing is diffuse and the mean free path is constant, such 
that the phonon conductivity hip oc T 3 . (Phonons can 
also be scattered by twin boundaries and grain bound- 
aries.) If the sample boundaries are smooth, specular re- 
flection yields a temperature-dependen t mea n free path, 
and n p (xT a , typically with 2 < a < 3pESl 

In Fig. [5j we show the thermal conductivity k c of our 
c-axis samples, plotted as k c /T vs T 2 , for magnetic fields 
from H = 0.05 to 15 T. Below T ~ 0.15 K, the curves 
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FIG. 5: Temperature dependence of the c-axis thermal 
conductivity « C) plotted as k c /T vs T 2 , for six samples of 
Ba(Fei_ I Co a: )2As2, with x as indicated, in a magnetic field 
H = 0.05, 4, and 15 T (data taken at other fields are not 
shown for clarity). These are six of the seven samples la- 
belled A in Table I. The lines are a linear fit to the data 
below T 2 — 0.015 K 2 , used to extract the residual linear term 
k c o/T as the extrapolation of k c /T to T = 0. The values of 
Kco/T are listed in Table I for H — 0, and plotted vs H in 
Fig. [6j for all twelve c-axis samples. Solid black squares on 
the T — axis (dashed line) give the residual linear term in 
the normal-state thermal conductivity, Kc,n/T, obtained from 
the residual resistivity p c o of the sample via the Wiedemann- 
Franz law (see Table I). 



are linear, consistent with diffuse phonon scattering on 
the sample boundaries of our c-axis samples, which are 
indeed characterized by rough side surfaces. We obtain 
K e /T = kq/T by extrapolating k/T to T — using a 
linear fit below T 2 — 0.015 K 2 . The error bar on this 
extrapolation is approximately ± 0.5 /iW/K 2 cm, for all 
c-axis samples. The value of k c q/T thus obtained is plot- 
ted as a function of field H in Fig. |6j for all twelve c-axis 
samples. For five concentrations, we have a pair of crys- 
tals with nominally the same Co concentration. As can 
be seen, the two curves in each pair are in good agreement 
with each other, well within the uncertainty in the geo- 
metric factor. The zero-field values are listed in Table I. 
They range from n c0 /T < 1 fiW /K 2 cm at x = 0.048 and 
0.074 to k c0 /T ~ 6 pW/K 2 cm at x = 0.038 and 0.127. 

The normal-state residual linear term k-^/T was es- 
timated using the values of po through application of 
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FIG. 6: Field dependence of the residual linear term k c o/T 
in the c-axis thermal conductivity of our twelve c-axis sin- 
gle crystals of Co-Bal22, with x as indicated. Underdoped 
compositions are shown on the left, overdoped compositions 
on the right. For five concentrations, two different crystals 
with nominally the same x value, labelled A (circles) and B 
(triangles), were measured (see Table I). 



the Wiedemann- Franz law. The value of Kn/T is shown 
as a solid black square on the y axis of Fig. [5] For 
the most heavily overdoped samples, with x — 0.127, 
a magnetic field of 15 T is sufficient to reach the nor- 
mal state, where kq/T saturates to its normal-state value 
kn/T. This allows us to check the Wiedemann- Franz 
law. For sample A, k c0 /T = 16.0 ± 0.5 pW/K 2 cm at 
H = 15 T, while k n /T = 16.3 ± 0.8 pW/K 2 cm; for 
sample B, n c0 /T = 20.0 ± 0.5 pW /K 2 cm at H = 15 T, 
while kn/T = 21.7 ± 1.1 pW /K 2 cm. Within error bars, 
associated with extrapolations to get kq/T and po, the 
Wiedemann- Franz law is satisfied in both samples. 



B. Heat transport in the a direction 

In Fig. [7j we show the thermal conductivity n a for six 
of our nine a-axis samples, plotted as n a /T vs T, for 
magnetic fields from H = to 15 T. Unlike in the c- 
axis samples, the phonon conductivity n p does not obey 
Kp/T oc T 2 as T — > 0. Instead, it follows approximately 
a power law such that k p /T oc T q_1 , with 2.0 < a < 2.5. 
These values of a are typi cal o f specular reflection off 
smooth mirror-like surfaces! 42 * 45 ! The cleaved surfaces of 
these Co-Bal22 crystals (normal to the c axis) are in- 
deed mirror-like. Previous measurements of in-plane heat 
transport on K-Bal22j23 Co-Bal22j^ and Ni-Bal22^1 
have all obtained a < 2.7. 

As done previously for other a-axis samples}^ we ob- 
tain the residual linear term k q0 /T by fitting the data 
below T = 0.3 K to a power-law expression, k/T = a + 
6r a_1 , where a = K a o/T. The error bar on this extrapo- 
lation is approximately in the range ± 10—20 /iW/K 2 cm. 
(The uncertainty is an order of magnitude larger than 
for k c q/T because the phonon-related slope is an order 
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FIG. 8: Field dependence of the residual linear term n a o/T 
in the a-axis thermal conductivity of our nine a-axis single 
crystals of Co-Bal22, with x as indicated. Underdoped com- 
positions are shown of the left, overdoped compositions on the 
right. For three concentrations, two different crystals with 
nominally the same x value, labelled A (circles) and B (tri- 
angles), were measured (see Table II). 



normalized to their respective normal-state value kn/T. 



FIG. 7: Temperature dependence of the a-axis thermal 
conductivity K a , plotted as K a /T vs T, for six samples of 
Ba(Fei_ :E Co a ;)2As2, with x as indicated, in a magnetic field 
H = 0, 4, and 15 T (data taken at other fields are not shown 
for clarity). These are the six samples labelled A in Table II. 
The lines are a power-law fit to the data below T = 0.3 K, 
namely k/T — a + fcT a_1 . The fit is used to extrapolate K a /T 
to T = and thus obtain the residual linear term K a o/T. The 
power q is in the range from 2 to 2.5. The values of K a o/T are 
listed in Table II for H = 0, and plotted vs H in Fig. [8] for 
all nine a-axis samples. Solid black squares on the T = axis 
give the residual linear term in the normal-state thermal con- 
ductivity, Kbn/T, obtained from the residual resistivity p a o of 
the sample via the Wiedemann- Franz law (see Table II). 



of magnitude steeper.) As found previously over the con- 
centration range 0.048 < x < 0.1 14p2l we again find 
Kao/T ~ 0, within error bars, now over a wider range: 
0.042 < x < 0.127. This is consistent with a separate 
report that K a0 /T ~ in Co-Bal22 at x = 0.135{22l 

Upon application of a magnetic field, K a o/T increases, 
as displayed in Fig. [8] for all nine a-axis samples. For 
three concentrations, we have a pair of a-axis crystals 
with nominally the same Co concentration. As can be 
seen, the two curves in each pair are in good agreement 
with each other, within the ± 20 % uncertainty in the 
geometric factor and the error bar on the extrapolations. 



IV. DISCUSSION 

The results of our study are summarized in Fig. |9j 
where the kq/T values of all 21 samples are plotted vs x, 



A. Gap nodes 

1. Zero magnetic field 

Our central finding is the presence of a substantial 
residual linear term kq/T in the thermal conductivity 
of Co-Bal22 in zero field, for heat transport along the c 
axis. It implies the presence of nodes in the supercon- 
ducting gap, such that A (k) = for some wavevectors 
k on the Fermi surfacePSEsEHH] Because heat conduc- 
tion in a given direction is dominate d by quasiparticles 
with k vectors along that direction) 34 ! 35 ! the fact that 
kq/T is negligible when heat transport is along the a 
axis, at all x, implies that the nodes are located in re- 
gions of the Fermi surface that contribute strongly to 
c-axis conduction but very little to in-plane conduction. 
The anisotropy of kq /T becomes pronounced as x moves 
away from the critical doping x s ~ 0.06, in either direc- 
tion. For x > 0.06, we see that the a/c anisotropy in 
ko/kn is at least a factor 10 (see Fig. [9]). Such a large 
anisotropy is not expected in a scenario of isotropic pair- 
breakingpS' and it confirms that the residual linear term 
seen in the c direction is due to nodes. 

At the highest doping studied here, x = 0.127, 
kq/ktss = 0.34 ± 0.03 for J || c. (This is for sample 
A, which has the lowest T c in the overdoped regime; 
see Table I.) This magnitude is typical of superconduc- 
tors with a line of nodes in the gap. In the heavy- 
fermion superconductor C elrlns , with T c — 0.4 K and 
H c2 ~ 0.5 T, k /k n ~ 0.2.^21421 i n the ruthenate super- 
conductor Sr 2 Ru04, with T c = 1.5 K and H c2 = 1.5 T, 
kq/^n — 0.1 — 0.3 (depending on sample purity)! 3 ^ In 
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FIG. 9: Residual linear term kq/T of Co-Bal22 normalized 
by the normal-state value kn/T as a function of Co concen- 
tration x, at H — (lower panel) and at H — i/ C 2/4 (upper 
panel). Full blue symbols are for a heat current along the c 
axis (J || c; circles for A samples, triangles for B samples in 
Table I). Empty red symbols are for a heat current along the 
a axis (J || a; circles for A samples, triangles for B samples 
in Table II) . The white interval between the two vertical grey 
bands at x < 0.032 and x > 0.17 is the region of supercon- 
ductivity in the phase diagram (Fig. [TJ. The vertical dashed 
line at x = 0.06 marks the approximate location of the critical 
doping x 3 where the structure at T = goes from orthorhom- 
bic (below) to tetragonal (above (see Fig.[l|. Lines through 
the data points are a guide to the eye. Error bars on the 
H = data are shown for the A samples (circles). 



the overdoped cuprate Ti2Ba2Cu06-<5 (Tl-2201), a d- 
wave superconductor with T c = 15 K and H c2 ~ 7 T, 
ko/^n — 0. 35P^1 In the latter case, because the or- 
der parameter is well-known and the Fermi surface is 
very simple (a single 2D cylinder), it was possible to 
show that the magnitude of k /T agrees quantitatively 
with the theoretical BCS expression for the residual lin- 
ear term in a d- wave sup erconductor,^ namely kq /T = 

where d is the interlayer sep- 
aration, fcp and Dp are the Fermi wavevector and ve- 
locity at the node, respectively, and S = SA/Sk is the 
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FIG. 10: Comparison of heat transport by nodal quasipar- 
ticles and the jump in heat capacity at the superconducting 
transition, in Co-Bal22 as a function of Co concentration x. 
Heat transport is measured as the zero-field residual linear 
term in the thermal conductivity along the c axis, k c o/T, 
normalized by the corresponding normal-state conductivity 
at T — > 0, k; c n/T, multiplied by T c . The heat capacity jump 
AC is divided by T c (from Ref. I5ip . The vertical dashed line 
marks the location of x s . Other lines are a guide to the eye. 



slope of the gap at the node. (For a d-wave gap with 
A(fc) = A cos(2<^), S = 2A .) 

If the line of nodes in the gap is imposed by the sym- 
metry of the order parameter, as in a d-wave state, then 
kq/T is universal, i.e. independent of the imp urity scat- 
tering rate Tq, in the clean limit HTq << Ao' 35 * 36 ^ Such 
universal transport was demonstrated experimentally for 
CeMn s Sa Sr 2 Ru0 4 , 1151 and the cuprates YBaaCuaO^ 
and Bi2Sr2CaCu2 08.^ As a fraction of the normal-state 
conductivity, one then gets (kq/T)/(k-^/T) = ko/ k n oc 
HTq/S^ However, if the nodes are not imposed by sym- 
metry, but are 'accidental', as in an 'extended-s-wave' 
state, they still cause a non-zero residual linear term, 
with Kq/T oc 1/S, but Kq/T is no longer universal, be- 
cause S depends on the scattering rate ToP^ 

In Fig. [9] we see that k c0 /k c n exhibits a strik- 
ing U-shaped dependence on Co concentration x, with 
Kco/^cN — > as x — > x s . Just above x s , at x = 0.074, 
k c0 /T = 0.2 ± 0.5 fiW/K 2 cm. (This is for sample B, 
which is closest to x s , as it has the highest T c ; see Ta- 
ble I.) This is equal to zero within error bars, indicating 
that there are no nodes in the gap at this concentration, 
as also inferred from the field dependence (see below) . If 
the nodes can be removed simply by changing x, then 
these nodes must be accidental, not imposed by symme- 
try. 

Given that the change in ko/kn with x on the over- 
doped side is due to a change in kq/T and not a change 
in kn/T (since po is independent of x, within error bars), 
we attribute the dramatic rise in «o/ k n from x — 0.074 
to x = 0.127 to a decrease of the slope S with increas- 
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ing x. Part of this decrease must be due to a drop in 
the overall strength of superconductivity, as measured 
by the decreasing T c . We can factor out that effect by 
multiplying kq/k-n by T c , as shown in Fig. 10 We see 
that aco/kn xT c vs x is far from constant, as it would be if 
the decrease of A(k) vs x was uniform, independent of k. 
In a (i-wave superconductor, for example, S would typi- 
cally scale with the gap maximum Ao, which itself would 
scale with T c , giving a constant product kq/k^x T c (for a 
constant Tq). By contrast, in Co-Bal22 the slope of the 
gap at the nodes decreases faster than that part of the 
gap structure which controls T c . In other words, A(k) 
must be acquiring a stronger and stronger k dependence, 
or modulation, with increasing x. 

In the underdoped regime, for samples with x — 0.048 
and lower, the metal is antiferromagnetic^ and its Fermi 
surface is reconstructed by the antiferromagnetic order. 
Nevertheless, a residual linear term k c q/T is still observed 
at H — (see Fig.|9|. At x — 0.038, it is even larger than 
at x = 0.127, namely k q /k n = 0.48 ± 0.04 (see Table I). 
This implies that nodes are present in the superconduct- 
ing gap inside the region of co-existing antiferromagnetic 
order. The fact that kq/T is again strongly anisotropic 
(see Fig. [9]) means that those nodes are still located in 
regions of the Fermi surface that contribute strongly to o 
axis conduction and little to a-axis conduction. The fact 
that the nodes survive the Fermi-surface reconstruction is 
consistent with their location in regions with strong 3D 
character, since the spin-density wave gaps the nested 
portions of the Fermi surface, which are typically those 
with strong 2D character. (It should be emphasized that 
the mechanisms responsible for the drop in T c and the 
rise in kq/k^ are likely to be different above and below 
optimal doping.) 



2. Field dependence 

The effect of a magnetic field on n n / T reveal s how easy 
it is to excite quasiparticles at T = o! 39 l 47 l 53 l For a gap 
with nodes, the rise in kq/T with H is very fast, because 
delocalized quasiparticles exist outside the vorticesp^ as 
shown for the d-wave superconductor Tl-2201 in Fig. [TT] 
For a full gap without nodes or deep minima, such as 
in the s-wave superconductor Nb, the rise in kq/T vs H 
is exponentially slow 



(see Fig. 11), because it relies on 
tunneling between quasiparticle states localized on ad- 
jacent vortices. For Co-Bal22 at x = 0.127, k c o/«n is 
seen to track the ci-wave data all the way from H = to 
H = H C 2 . This nicely confirms the presence of nodes in 
the gap structure of overdoped Co-Bal22 that dominate 
the transport along the c axis. 

By contrast, at x = 0.074, the initial rise in k c q/T vs 
H has the positive (upwards) curvature typical of a node- 
less gap, for both samples A and B (see Fig. [6]). The rise 
at low H is faster than in a simple s-wave superconduc- 
tor like Nb (Fig. Ill, either because of a fc-dependence 
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FIG. 11: Bottom panel: Residual linear term Ko/T of Co- 
Bal22 normalized by the normal-state value kn /T as a func- 
tion of magnetic field H, plotted as (ko/T)/(kn/T) = ko/kn 
vs H/H C 2 for three representative Co concentrations, as in- 
dicated: underdoped (x = 0.042; red), slightly overdoped 
(x = 0.074; black), and strongly overdoped (x = 0.127; blue). 
kn/T is obtained from the Wiedemann- Franz law (see text 
and Tables I and II), except for the samples with x = 0.127 
where we use the value of kq/T measured at H = 15 T, since 
H C 2 — 15 T at that concentration. Full circles are for a heat 
current along the c axis (J || c; data from three A samples - 
see Table I). Empty circles are for a heat current along the a 
axis (J || a; data from three A samples - see Table II). The 
vertical dashed line marks H = _H" C 2/4; the value of ko/^n at 
i/ C 2/4 is plotted in the top panel of Fig. [9] for all x. We also 
reproduce corresponding data for the d-wave superconductor 
Tl-2201 (from Ref. [50}, the isotropic s-wave superconductor 
Nb, and the multi-band s-wave superconductor NbSe2 (from 
Ref. I54[) . Top panel: Anisotropy of the normalized residual 
linear term ko/kn at x = 0.127. The red dashed line at 
H C 2/10 marks roughly the field beyond which ko/ k n becomes 
isotropic. 



of the gap or because of a multi-band variation of the 



gap amplitude, or both. A multi-band variation is what 
causes the fast initial rise in kq/T vs H (with positive 
curvature) in NbSeJ^l ( se e Fig. 11 ). This H dependence 
strongly suggests that there are no nodes in the gap of 
Co-Bal22 at x = 0.074, as inferred above from the neg- 
ligible value of K a o /T. 
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B. Gap minima 

We saw that nodes in the gap have two general and re- 
lated signatures in the thermal conductivity 1) a finite 
residual linear term Kq/T in zero field, and 2) a fast initial 
rise in kq/T with H. Both signatures are clearly observed 
in Co-Bal22 at x = 0.127 for J || c. For J || a, however, 
the situation is quite different. Indeed, n a o/T is negligi- 
ble at H = 0, for all x, as also found in previous measure- 
ments of K a on underdoped K-Bal22, 21 optimally-doped 
Ni-Bal22,S£l and overdoped Co-Bal22P 

Consequently, the fast initial rise in kq/T with H for 
J || a, seen in Fig. is not due to nodes but rather 
to the presence of deep minima in the gap, in regions 
of the Fermi surface that contribute significantly to in- 
plane conduction, as previously reported.^! In the top 
panel of Fig. [9j we show the normalized residual linear 
term ko/ k n measured at H = -ff C 2/4. We see that in 
the overdoped regime the ko/ k n values are the same 
for both current directions, at all x. In other words, 
whereas «o/ k n is very anisotropic at H = 0, it is essen- 
tially isotropic at H > 7J C 2/10, as shown for x = 0.127 in 



the top panel of Fig. 11 But quasiparticle transport for 
J 1 1 c is due to nodal excitations, whereas quasiparticle 
transport for J || a comes from field- induced excitations 
across a minimum gap. In a single-band model, say with 
a single ellipsoidal Fermi surface, this contrast between 
zero-field anisotropy and finite-field isotropy can only be 
described by invoking two unrelated features in the gap 
structure A(k): nodes along the c axis and deep minima 
in the basal plane. However, the fact that ko/kn remains 
isotropic at all x (for H — H c2 /A) strongly suggests that 
nodes and minima are in fact intimately related. We 
therefore propose that they both come from the same 
tendency of the gap function A(k) to develop a strong 
modulation as a function of k, which causes a deep mini- 
mum on one Fermi surface and an even deeper minimum 
on another Fermi surface, where the gap would actually 
go to (or through) zero. In other words, instead of in- 
voking two unrelated features of the gap structure on a 
single Fermi surface, we invoke a single property of the 
gap structure which leads to two related manifestations 
on separate Fermi surfaces. 



C. Two simple models for the gap structure 

For the purpose of illustration, we consider a simplified 
two-band model for the Fermi surface, whereby one sur- 
face has strong 3D character and the other has quasi-2D 



character, as sketched in Fig. 12 The 3D Fermi surface 
can either be open along the c axis, as drawn in Fig. [12] 
and suggested by some ARPES dataPor closed, as sug- 
gested by some band structure calculations.^ The 3D 
Fermi surface is responsible for most of the c-axis conduc- 
tion and the 2D surface for most of the a-axis conduction 
(recall that in the T = limit K a N/K C N = Pcn/p&n — 20). 
Note that in reality the Fermi surface of Co-Bal22 con- 



tains at least four separate sheets EMU our model requires 
that at least one of these has strong 3D character and it 
treats all others in terms of a single Fermi surface, the 
second quasi-2D sheet. We then propose that the gap 
A(k) varies strongly as a function of k, on both Fermi 
surface sheets. There are two basic scenarios: a gap mod- 
ulation as a function of k z , illustrated in Fig. |12[ or a gap 
modulation as a function of the azimuthal angle <p in the 
basal plane, illustrated in Fig. 13 The strong modula- 



tion extends to negative values on the 3D Fermi surface, 
thereby producing nodes where A(k) ~ 0, whereas it 
only produces a deep minimum (where A = A m i n ) on 
the 2D Fermi surface (at least in the range of concentra- 
tions covered here). In the first scenario (Fig. 12 1, the 
lines of nodes are horizontal circular loops in a plane nor- 



mal to the c axis; in the second scenario (Fig. 13), they 
are vertical lines along the c axis. 

Both versions of the model explain the isotropy at 
H C 2/A and the anisotropy at H = 0. The isotropy of 
ko/^n follows fundamentally from having a similar k 
modulation of the gap on both Fermi surfaces. When 
the field is large enough to excite quasiparticles across 
the minimum gap on the 2D Fermi surface, quasiparti- 
cle transport from both Fermi surfaces will be similar, 
explaining the rapid and isotropic rise in kq/T with H. 
By contrast, at H = no quasiparticles are excited on 
the 2D Fermi surface at T = (since fceT << A m j n ), 
whereas nodal quasiparticles are always present on the 
3D surface. This explains the large anisotropy of ko/kn 
at H = 0. Note that this anisotropy is not governed 
by the anisotropy of the gap itself, i.e. by the direction 
of the gap modulation, but rather by the fact that the 
nodes lie on the 3D Fermi surface. (Whether horizontal 
or vertical line nodes are more consistent with our data 
depends on details of the real Fermi surface of Co-Bal22.) 
The nodal quasiparticles on the 3D sheet must also con- 
tribute to a-axis conduction. Assuming that for the 3D 
Fermi surface K a o/T — k c q/T at H — 0, we should detect 
a residual linear term n a o/T ~ 6 /iW/K 2 cm in the a- 
axis sample with x — 0.127, for example. This is indeed 
consistent, within error bars, with the value we extrap- 
olate for the a-axis data at x = 0.127 (Fig. |7j), namely 
K a0 /T = 17 ± 20 ^W/K 2 cm (Table II). 

In both versions of our model for the gap structure, 
the U-shaped x dependence of Kq/k^ is attributed to an 
increase in the modulation of the gap as x moves away 
from x s , as illustrated in Figs. [12] and 13 The fact that 
the U-shaped curves in Fig. [9] have their minimum where 
the (inverted U-shaped) T c vs x curve has its maximum 
points to a reverse correlation between T c and gap mod- 
ulation. Modulation is a sign of weakness. The presence 
of nodes in the gap may then be an indicator that pairing 
conditions are less than optimal. 

It is possible that at high enough x in the overdoped 
regime A m j n , the minimum value of the gap on the quasi- 
2D Fermi surface, goes to zero, so that nodes appear on 
that Fermi surface as well. This would immediately cause 
Kao/T to become sizable. It is conceivable that the large 
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FIG. 12: First model for the gap structure. Simplified two- 
band model of the Fermi surface of Co-Bal22 (blue solid line), 
shown in the a — c plane, with k z j c in the vertical direction. 
Although in reality the Fermi surface of Co-Bal22 consists of 
at least four sheets, in our model we reduce it to two sheets: a 
sheet with strong 3D character (FS # 1; top left) and a sheet 
with quasi-2D character (FS #2; top right). The supercon- 
ducting gap A(k) on both Fermi surfaces (red dashed line) 
varies strongly as a function of k z , as shown on the bottom 
for three representative concentrations x from x s upwards. 
On the 3D Fermi surface (#1), the gap modulation is such 
that it extends to negative values, producing nodes (black 
circles) at certain points. On the 2D Fermi surface (#2), the 
gap modulation is strong enough to cause a deep minimum 
(where A = A m i n ), but not nodes. At x — x 3 , the gap minima 
are shallow and there are no nodes on either Fermi surface. 
With increasing x (decreasing T c ), the modulation increases, 
the minima deepen and the nodes appear. A further increase 
in x (beyond the maximal concentration in this study) could 
eventually also yield nodes on FS #2. 



value of K a Q/T measured in undoped KFe2As2,^ which 
can be viewed as the strong doping limit of K-Bal22, is 
the result of a gap modulation so strong that it goes to 
(through) zero on all Fermi surfaces. 

A pronounced modulation of A(k) should manifest it- 
self in a number of physical properties. For example, in 
an s-wave superconductor, a variation of the gap magni- 
tude over the Fermi surface, whether from band to band 
as in MgB2p^or from k dependence (anisotropy) as in 
ZnJ^I leads to a suppressed ratio of specific heat jump 
AC at the transition to T c . The pronounced gap mod- 
ulation and anisotropy revealed by the thermal conduc- 
tivity could therefore account for the dramatic variation 
of AC/T C measured in Co-Bal22 vs reproduced in 
Fig. 10 AC/T C is seen to be maximal where kq/kn is 



FIG. 13: Second model for the gap structure. Same as in 
Fig. |12| except that the modulation of the superconducting 
gap A(k) is now a function of the azimuthal angle (f> in the 
basal plane. 



minimal, i.e. where the gap modulation is weakest, and 
it drops just as rapidly with a change in x as kq/^n rises. 



D. Theoretical calculations 

The two-band picture suggested by our thermal con- 
ductivity data is reminiscent of the proximity scenario 
proposed for the three-band quasi-2D p-wave supercon- 
ductor Sr 2 Ru04f^l where superconductivity originates 
on one band, the most 2D one, and is induced by prox- 
imity on the other two bands. This fc-space promixity 
effect is such that a k z modulation of the induced gap 
produces horizontal line nodes on the latter two Fermi 
surfaces,^! in analogy with the horizontal-line scenario of 
Fig. 12 A proximity scenario of this sort was in fact pro- 
posed for the pnictidesp^ predicting c-axis nodes in the 
superconducting gap. The effect on the superconducting 
gap structure of including the k z dispersion of the Fermi 
surface in BaFe2As2 was recently calculated within a 



spin-fluctuation pairing mechanism on a 3D multi-orbital 
Fermi surface.^! A strong modulation of the gap A(k) as 
a function of both k z and <f> is obtained which can indeed, 
for some parameters, lead to accidental nodes. 

The thermal conductivity of pnictides was calculated 
in a 2D two-band model for the case of an extended-s- 
wave gap (of A\ g symmetry) P3 These calculations show 
that the presence of deep minima in the gap, in this case 
as a function of (j), can account for the rapid initial rise 
observed in n a o/T v s H, starting from K a o/T — at 
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H = 0. It seems clear that calculations for a gap whose 
deep minima occur instead as a function of k z would yield 
similar results. It will be interesting to see what calcu- 
lations of the thermal conductivity give when applied to 
the 3D model of Ref . or indeed to the simple two-band 
models proposed here (in Figs. 12 and 13). 



V. CONCLUSIONS 

In summary our measurements of the thermal 
conductivity in the iron-arsenide superconductor 
Ba(Fei_ x Coj;)2As2 show unambiguously that the gap 
A(k) has nodes. These nodes are present in both the 
overdoped and the underdoped regions of the phase 
diagram, implying that they survive the Fermi-surface 
reconstruction provoked by the antiferromagnetic order 
in the underdoped region. The nodes are located in 
regions of the Fermi surface that dominate c-axis conduc- 
tion and contribute very little to in-plane conductivity. 
The fact that the strongly anisotropic quasiparticle 
transport at H — becomes isotropic in a magnetic field 
H = H c2 /4: shows that there must be a deep minimum 
in the gap in regions of the Fermi surface that dominate 
in-plane transport. These two features - nodes on 3D 
regions and minima on 2D regions of the Fermi surface 
- point to a strong modulation of the gap as a function 
of k. This modulation of A(k) would be present on all 
Fermi surfaces, but be most pronounced on that surface 
with strongest k z dispersion, where it has nodes. This 
suggests a close relation between the 3D character of the 
Fermi surface and gap modulation. 

The anisotropy of k shows a strong evolution with Co 
concentration x. At optimal doping, where T c is maxi- 
mal, there are no nodes and kq/T has the anisotropy of 
the normal state. With increasing x, nodes appear and 
kq/T acquires a strong anisotropy. We attribute this to 
an increase in the gap modulation with x, which may 
explain the strong decrease in the specific heat jump at 
TjPAl and the change in the power-law temperature depen- 



dence of the penetration depthl^E^ The fact that nodes 
are located in regions that dominate c-axis conduction is 
consistent with the fact that the penetration depth along 
the c axis has a linear temperature dependence.^ 

Horizontal line nodes in Co-Bal22, which would be 
the result of a strong modulation of the gap along k z 
rather than a strong in-plane angular dependence, would 
reconcile the isotropic azimuthal angular dependence of 
the gap seen by ARPES with the evidence of nodes or 
minima from thermal conductivity, NMR relaxation rate 
and penetration depth measurements in Co-Bal22 and 
other iron-based superconductors. A k z modulation of 
A(k) should be detectable by ARPES, especially in the 
overdoped regime where it would be strongest. 

Because the nodes go away by tuning x towards opti- 
mal doping, we infer that they are 'accidental', i.e. not 
imposed by symmetry, and so consistent a priori with 
any superconducting order parameter, including the s± 
stateP^HEl Although accidental nodes are not a direct 
signature of the symmetry, the strong modulation of the 
gap nevertheless reflects an underlying k dependence of 
the pairing interaction, and as such the 3D character of 
the gap function A(k) is an important element in under- 
standing what controls T c in this family of superconduc- 
tors. 
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